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INTRODUCTION 


The  time  and  frequency  domains  are  equivalent  displays  of  seismic  trace  information,  though 
some  qualities  of  the  signal  are  more  easily  observed  in  one  domain  than  the  other.  The  relative 
frequency  excitation  of  Lg,  for  instance,  is  most  easily  viewed  in  the  frequency  domain,  but  such 
waveform  qualities  as  the  sequence  in  which  pulses  arrive  in  the  wave  train  or  the  sharpness  of 
pulse  onset  are  most  easily  studied  in  the  time  domain  (Murphy  and  Bennett,  1982;  Blandford, 
1981).  Because  of  the  tremendous  complexity  of  high  frequency  regional  data,  most  attempts  at 
using  it  for  discrimination  purposes  have  involved  analysis  of  the  frequency  content  of  the  various 
arrivals  either  through  transforming  selected  windows  or  through  multiple  bandpass  filtering.  We 
report  here  on  our  initial  attempts  to  explore  the  alternative  and  to  discriminate  events  using  those 
waveform  characteristics  most  easily  observed  in  the  time  domain. 

A  second  advantage  of  time  domain  analysis  approaches  is  that  they  permit  a  deeper  insight 
into  the  physical  processes  creating  a  seismic  signal's  character.  For  this  reason,  they  can  be  more 
easily  used  to  evaluate  the  transportabilty  of  a  discriminant  to  varying  geophysical  and  tectonic 
regimes.  This  is  an  especially  important  feature  in  the  development  of  regional  discriminants.  The 
most  prominent  and  successful  spectral  regional  discriminants  have  been  empirically  developed. 
This  means  that  they  must  be  redeveloped  and  reverified  ir.  each  new  area.  As  we  shall  show  in 
the  lollowing,  through  rigorous  time  domain  analysis  such  features  as  regional  depth  phases  can 
be  identified  and  used  to  discriminate.  Discriminants  based  on  such  simple  physical  features  as 
source  depth  should  be  transportable  anywhere. 

In  work  recently  completed  under  the  treaty  verification  program,  we  have  proved  that  such 
time  domain  discriminants  do  exist.  In  analyzing  a  test  discrimination  data  set  from  the  western 
U.  S.,  we  have  discovered  that  the  onset  of  Pn  is  always  very  similar  for  explosions  and  that  few 
earthquakes  have  this  unique  waveform  character.  This  information  can  be  constructed  into  a 
simple  discrimination  scheme  by  testing  the  correlation  of  observed  Pn  waveform  onsets  with 
average  waveforms  observed  from  explosions.  High  correlations  indicate  explosions  and  low 
correlations  earthquakes.  We  have  also  discovered  that  the  regional  phase  Pg  is  actually  composed 


1 


of  a  sequence  of  sub-arrivals  which  correspond  to  successively  higher  orders  of  reverberation  in 
the  crust.  In  realistic  crust  models,  the  depth  phases  play  an  important  rol-'  in  the  waveshapes  of 
these  sub-arrivals.  By  selecting  an  appropriate  frequency  band  to  analyze,  we  have  been  able  to 
accurately  model  this  type  of  data  from  explosions  in  the  western  United  States.  Oser  the  very 
relevant  regional  distance  ranges  of  200  to  600  km,  it  appears  that  a  discrimination  procedure  very 
similar  to  the  one  which  is  known  to  work  for  Pn  will  also  be  effective  for  Pg.  We  are  investigating 
whether  similar  discriminants  can  be  constructed  based  on  the  phases  Sn  and  Sg  in  areas  where 
those  phases  are  prominent  arrivals. 


RECENTLY  DEV  FLO  PEP  WAVEFORM  DISCRIMINANTS 
Though  the  technology  for  recording  broad  band  seismic  data  digitally  has  existed  for  some 
time,  a  good  regional  net  of  stations  surrounding  NTS  has  only  been  put  into  place  recently  (Figure 
I).  Several  different  types  of  stations  and  seismometers  are  in  the  net  including  DWWSSN  (JAS. 
AI.Q  or  ANMO),  LLNL  broad  bands  (EL.K,  MNV,  LAC',  KNHiand  the  new  Streckeisens  installed 
by  Caltech  (PAS)  and  UCSD  (PFO).  Much  recent  effort  at  Woodward  Clyde  has  been  directed  at 
developing  short  period  regional  discriminants  that  work  on  data  from  this  net.  The  situation  of 
having  data  available  from  many  sources  veiy  close  to  each  other  which  are  similar  in  character 
such  as  NTS  explosions  is  a  unique  one.  Fhe  timing  and  locations  of  the  events  are  known  exactly. 
This  turns  any  single  regional  station  in  the  western  U.  S.  into  the  equivalent  of  a  regional  seismic 
array.  Any  of  the  standard  array  processing  techniques  can  be  used  with  the  role  of  sources  and 
stations  being  reversed.  There  is  some  variability  in  source  time  historv  and  neat  source  structure, 
but  on  the  other  hand  the  receiver  structure  is  constant.  We  began  our  previous  work  by  simply 
obtaining  suites  of  records  from  each  of  the  stations  in  the  net  and  summing  to  form  stable  averages. 
In  so  doing,  we  observed  important  features  in  the  average  Pn  and  the  average  Pg  waveforms  which 
could  each  be  used  to  develop  time  domain  discriminants.  Our  success  to  date  is  summarized 
below. 

The  Pn  Waveform  Discriminant  Our  initial  success  in  developing  a  Pn  waveform  discriminant  was 
rooted  in  past  experience  with  modeling  waveforms  of  teleseismic  P  waves  from  nuclear  explosions. 
Figure  2  is  taken  Fom  Burdick  et  al.  (1984).  It  shows  observed  and  synthetic,  long  and  short 
period  records  from  the  nuclear  test,  CANNIKIN.  Arrows  on  the  left  of  the  figure  draw  attention 
to  a  subtle  feature  in  the  short  period  records  which  is  associated  with  the  arrival  of  the  phase  pP. 
The  records  where  the  feature  appears  are  assumed  to  be  along  high  Q  paths  It*  =  0.8  s).  Along 
low_*r  Q  paths  the  feature  washes  out.  The  long  periods  are  not  affected  by  pP  in  a  clearly  visible 
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Figure  2.  Comparison  ol  observed  short  and  long  period  P  waves  with  svntheuc  seismograms  for 
a  range  of  t  values  at  teleseismic  distances.  As  t*  increases,  the  interference  in  the  second  upswing 
produced  by  pP  becomes  less  apparent.  pP  can  be  observed  as  a  distinctive  double  peak  and  is 
produced  in  the  synthetics  with  an  artificial  delav. 


This  previous  experience  with  the  effects  of  pP  on  teleseismic  short  period  P  waves  proved 
valuable  in  interpreting  regional  P„  waveforms.  As  it  happens,  some  of  the  stations  in  the  western 
U.S.  digital  net  have  a  band-limited  response  equivalent  to  the  WWSSN  stations  in  Figure  2.  To 
make  a  uniform  comparison,  all  of  the  records  in  the  data  base  were  translormed  to  short  perioc 
WV.'SSN  instrument  response.  When  the  digital  signals  from  NTS  explosions  were  averaged,  it  was 
observed  that  the  average  Pn  waveform  was  very  similar  to  those  on  the  left  of  Figute  2  although 
much  shorter  in  period.  Average  waveforms  at  5  stations  in  the  digital  net  are  shown  in  Figure 
3  The  consistent  splitting  of  th^  second  upswing  is  indicated  by  the  arrows.  The  cleat  implication 
is  that  the  physics  of  the  wave  piopagation  of  short  period  teleseismic  P  and  short  period  Pn  is 
very  comparable.  More  precisely,  the  interaction  of  pP  with  P  must  be  similar  in  the  two  instances. 
This  is  not  unreasonable  in  that  the  apparent  velocity  of  Pn  is  about  8  km,  s  and  that  ot  teleseismic 
P  only  incteases  to  about  12  kims  at  30°  The  associated  change  in  pP  timing  and  amplitude  is 
small  The  implications  for  regional  discrimination  are  clear.  Only  very  shallow  sources  like 
explosions  will  have  depth  phases  at  very  short  times  Fatrhquake  depth  phases  will  be  much  later 
To  test  the  performance  of  this  discriminant,  we  assembled  a  set  of  Pn  waveforms  front  small 
earthquakes  near  NTS,  windowed  out  the  first  three  seconds  of  Pn  and  measured  the  correlation 
with  the  average  Pr,  waveforms  of  explosions  like  those  shown  in  Figure  3  A  similar  procedure 
was  carried  out  on  the  explosion  data  base.  The  results  from  one  typical  station,  J.AS  are  shown 
in  Figure  4  The  explosions  are  displayed  as  stars  and  the  earthquakes  as  crosses  The  separation 
of  th»  pope'nons  is  good  enough  to  warrant  more  study  of  this  discriminant. 

1  he  discrimination  capability  illustrated  in  Figure  4  only  demonstrates  that  the  Pn  waveform 
of  explosions  is  stable  at  JAS  and  consistently  different  from  earthquake  waveforms.  However, 
the  similarity  of  the  waveforms  in  Figure  3  suggests  that  the  shape  of  the  explosion  waveform  is 
consistent  from  station  to  station.  That  this  is  indeed  the  case  is  shown  in  Figure  5  where  the 
average  JAS  explosion  waveform  has  been  correlated  with  the  explosion  data  base  at  MNV. 


Figure  3.  The  average  Pn  waveform  for  Pahutc  explosions  observed  at  stations  in  the 
U.S  digital  net.  The  split  third  swing  is  indicated  by  arrows. 
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i-igurc  5.^  Discrimination  ol  Yucca  explosion:;  I'rom  earthquakes  using  the  average  Pn  waveform 
‘ rorp  and  .he  data  base  from  JAS.  Explosions  arc  stars  and  earthquake  are  crosses.  The 

separation  of  the  populations  is  oomparacle  to  that  m  Figure  4  illustrating  the  transportability,  of 
the  discriminant. 


Moderately  good  event  discrimination  is  still  possible.  These  studies  suggest  that  one  method  to 
achieve  discrimination  with  a  regional  net  will  be  to  continuously  monitor  the  waveform  of  Pn 
onset  at  all  of  the  stations  and  to  test  in  an  objective  fashion  for  the  presence  of  depth  phases. 

The  small  shoulder  on  the  second  upswing  of  the  waveforms  in  Figure  3  seems  to  be  :■ 
relatively  subtle  feature  on  which  to  base  a  discriminant.  However,  its  apparently  subile  nature 
is  due  to  the  particular  response  of  the  WW'SSN  instrument  To  demonstrate  this,  we  carried  out 
a  forward  modeling  study  of  the  broad  band  Pn  waveforms  from  the  LLNL  stations  designed  to 
determine  the  true  size  of  the  pP  arrival.  A  typical  result  is  shown  in  Figure  6  where  the  waveforms 
of  four  Pahute  events  as  recorded  at  MNV  are  displayed.  The  data  are  shown  as  solid  lines  and 
synthetics  for  an  appropriate  crustal  model  and  explosion  source  are  shown  as  dashed  lines.  Because 
the  instruments  are  broad  band,  deconvolving  the  response  out  is  a  stable  operation  which  has  been 
carried  out.  The  traces  shown  are  true  ground  velocity.  The  arrows  indicate  the  arrival  of  pP., 
in  the  data  and  synthetics.  The  observed  pPr  is  consistently  much  later  than  the  elastic  predictions. 
This  discrepancy  was  observed  for  all  events  at  all  stations  where  the  waveforms  could  be  modeled. 
In  many  instances,  it  was  observed  that  pP  was  as  much  as  two  times  larger  than  the  elastic 
predictions  Similar  results  (pP  late  and  amplified)  have  been  reported  in  most  studies  where  pP 
times  from  .  uclear  tests  have  been  measured.  In  our  final  report  for  our  previous  work  (Burdick 
et  al  .  1988),  we  showed  that  this  phenomenon  could  be  explained  in  terms  of  spall.  These  studies 
are  relevant  here  in  two  ways.  First,  effective  pP  ^that  is  pP  plus  spall)  is  a  large  arrival,  and  u 
is  sensible  to  attempt  to  continuously  examine  signal  onsets  for  its  presence  in  order  to  discriminate 
Second,  synthetic  modeling  of  the  P„-pPM  interaction  is  very  feasible,  and  it  is  reasonable  to  attempt 
its  use  in  automated  discrimination  schemes. 

The  Pg  Discriminant  The  c  ;.id  new  short  period  regional  discriminant  we  are  developing  and 
testing  in  the  western  L’  based  on  the  observation  that  the  Pg  phase  turns  out  to  contain  a 

sequence  of  sub-arrivals.  They  coirespond  to  successive  reverberations  of  energy  in  the  crustal 
wave  guide,  and  we  have  named  them  crustal  resonance  phases.  Our  initial  evidence  as  to  then 
existence  and  character  emerged  when  we  averaged  observations  of  suites  of  events  al  a  single 
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Figure  6.  Observed  versus  synthetic  Pn  waveforms  at  MNV  for  Pahute  events.  The  turning  ray 
model  for  propagation  is  assumed.  The  Mueller-Murphy  source  scaling  relations  and  a  i'  value  of 
0.1  s  are  also  assumed.  As  indicated  by  arrows,  the  observed  pP  arrival  is  iate  but  approximately 
the  same  in  size  as  the  elastic  prediction. 


station  just  as  did  the  evidence  for  pPn.  Though  there  were  some  indications  of  a  series  of  arrivals 
pumping  energy  into  the  Pg  coda  in  some  of  the  raw  data,  the  effects  of  scattering  and  othe; 
variations  between  events  seems  to  dominate.  This  has  long  been  known  to  be  the  greatest  difficulty 
in  interpreting  regional  data.  To  suppress  the  scattering  effects,  we  aligned  the  signals  very  carefully 
at  Pn  time  and  summed.  The  Yucca  and  Pahute  data  were  treated  separately.  This  procedure 
neglects  the  difference  in  slowness,  but  our  results  indicate  that  this  is  a  reasonable  approach 
Some  of  the  clearest  evidence  for  the  phases  came  from  ALQ.  The  res'Tt  for  Yucca  events  is 
shown  in  Figure  7  where  we  display  a  progressive  sum  of  records.  Each  subsequent  line  shows 
the  effect  of  adding  in  an  addition1.!  record  The  final  sum  clearly  indicates  that  there  are  coherent 
resonance  phases  at  frequencies  of  several  hz.  The  first  two  crustal  resonance  phases  develop 
almost  immediately  indicating  that  they  are  very  coherent.  A  third  appears  b'  the  time  4  traces 
are  added  in  and  the  sum  has  stabilized  by  the  time  8  are  added  m  Actually,  this  is  a  relatively 
rapid  rate  of  stabilization  indicating  that  the  resonance  phases  are  very  coherent  for  this  path.  At 
the  bottom,  we  show  a  synthetic  for  a  simple  layer  over  a  half  space  model  to  indicate  that  the 
resonance  phases  have  approximately  the  character  we  expect. 

Since  it  is  clear  that  these  reverberation  phases  exist,  it  is  natural  to  explore  the  possibility 
that  they  could  be  useful  in  discrimination.  There  have  been  many  studies  of  the  relative  properties 
of  the  long-duration,  composite  phases  such  as  Pg  or  Sg.  Now  that  we  know  that  the  complete 
phases  aie  built  up  of  sub  pulses,  we  can  attempt  to  base  discriminants  on  their  properties.  A 
large  number  of  possibilities  exist  Any  discriminant  which  has  been  tried  on  the  composite  pulses 
can  be  tested  on  the  sub  pulses.  These  could  include  particle  motion,  relative  amplitude,  frequency 
content  or  sharpness.  To  illustrate  that  important  variations  in  the  character  of  the  resonance 
phases  exist,  we  show  the  three  component  record  of  the  Pahute  event  flEKRA  at  ANMO  in 
Figure  8.  The  instrument  response  is  short  period  DWWbSN.  The  resonance  phases  are  clear  on 
the  vertical,  but  vary  markedly  iri  their  character  on  the  horizontals.  The  second  resonance  phase 
appears  to  he  carrying  much  more  S  energy.  As  we  will  show  in  the  following,  the  first  crustal 
resonance  appears  to  be  an  eff  icient  energy  channel  lor  P  and  later  resonances  for  S  One  possible 
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Figure  7.  The  crustal  resonance  phases  for  Yucca  events  are  separated  from 
the  noise  by  stacking.  The  synthetic  at  the  bottom  merely  indicates  that  the 
phases  have  approximately  the  right  timing  and  amplitude. 
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FLrure  3.  Three  component  observation  ac  AhXO  of  a  small  Pahute  event.  Noce 
the  increase  in  the  amount  of  S  energy  at  the  time  of  arrival  of  the  secon 
crustal  resonance  phase.  This  is  presumed  to  be  scattered  SV . 


method  for  discriminating  using  Pg  may  be  to  measure  the  relative  amount  cf  scattered  SV  energy 
between  the  first  and  third  resonance  phases.  There  should  be  more  scattered  S  in  earthquakes. 
In  the  following,  we  shall  discuss  the  development  of  a  time  domain  discriminant  xsed  on  time 
domain  correlations  of  longer  period  (2  sec.  and  greater)  energy,  but  the  point  here  is  that  there 
are  a  number  of  alternate  discriminants  possible. 

We  began  our  work  in  the  western  US  by  attempting  to  identify  a  frequency  band  in  which 
the  stable  resonance  phases  (resolved  in  Figure  7  by  stacking)  could  be  cbser\cd  consistently  in 
individual  records.  This  would  naturally  occur  in  a  somewhat  longer  period  band.  We  obtained 
suites  of  explosion  records  from  each  of  the  stations  in  the  digital  net  and  experimented  with  a 
number  of  bandpass  filters.  We  achieved  the  desired  stability  using  the  WWSSN  long  period  response 
modified  by  a  high  pass  and  a  low  pass  third  order  butterworth  filter.  The  long  period  butterworth 
cutoff  was  positioned  at  10  sec  and  the  low  pass  at  .6  sec.  This  resulted  in  signals  with  a  dominant 
period  of  about  2  sec.  A  suite  of  records  from  about  15  events  was  released  to  us  by  LLNL  and 
some  of  the  records  with  the  best  signal  to  noise  are  dispiayed  in  Figure  9.  The  strong  correlation 
of  the  waveforms  is  apparent  as  is  the  move  out  of  Pn  in  front  of  Pg  with  range.  There  are 
differences  between  the  Yucca  Flat  and  Pahute  Mesa  signals  at  the  same  station,  but  this  is  not 
surprising  given  the  variation  in  distance  to  the  test  sites 

We  developed  and  tested  appropriate  generalized  ray  and  wavenumber  integration  codes  for 
computing  synthetics  of  Pg  which  in  this  period  range  is  sometimes  called  Pnl.  We  examined  a 
number  of  plausible  crustal  models  attempting  to  keep  them  as  simple  and  realistic  as  possible 
Our  preferred  structure  model  for  the  western  U.S.  was  a  crustal  layer  32  km  thick  over  a  standard 
lid  model  with  a  realistic  free  surface  velocity  decrease  grading  smoothly  over  the  top  6  km.  Green's 
functions  for  the  model  computed  using  wavenumber  integration  and  generalized  ray  theory  are 
compared  on  the  left  and  right  of  Figure  10  respectively  for  a  suite  of  ranges.  The  wavenumber 
results  are  exact  except  for  the  limits  of  numerical  integration  while  the  generalized  ray  results 
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Figure  9.  Recorded  Pn j  seismograms  from  both  Yucca  Hat  and  Paiiute  Mesa  on  the  LLNL  network 
after  convolving  with  a  WWSSN  LP  instrument  and  a  recursive  bandpass  filter  i.low-cut=  0.1  Ha, 
high  cut  »  0.6  Ha  and  third-order  Butterworth  causal).  (,a)  For  Yucca  flat  explosions  and  tb)  tor 
Pahute  Mesa  explosions.  Data  show  remarkable  stability.  Original  instrument  response  was 
deconvolved. 
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Figure  10.  A  comparison  of  Green's  functions  computed  using  wavenumber 
integration  vs.  generalized  rays.  Traces  shown  are  displacement  with 
a  Mue 1 1 er -Murphy  source,  four  dominant  ray  groups  are  identified. 
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are  somewhat  affected  by  truncation  of  the  ray  sum.  The  generalized  rays  include  only  those 
which  reflect  from  the  free  surface  or  the  Moho  thus  neglecting  internal  reflections  within  the 
surface  gradient. 

The  utility  of  having  reasonably  accurate  generalized  ray  synthetics  is  that  they  can  be 
decomposed  into  ray  subsets  to  determine  the  paths  carrying  the  dominant  energy.  Four  subgroups 
of  rays  are  identified  in  the  figure,  and  they  are  clearly  related  to  the  resonance  phases  discussed 
above.  The  first  energy  is  of  course  Pn,  but  the  first  large  arrival  or  Pg  onset  is  composed  ot  PmP 
and  pPmP.  The  contribution  of  the  free  surface  reflection  is  significant  The  second  group  of 
rays  forming  the  second  resonance  is  dominated  by  2PmP  and  2pPmP  The  third  group  labeled 
C.  Gr  refers  to  some  unusual  rays  in  a  Converted  Group.  They  are  dominated  by  PmSPmP, 
PmPSmP  and  their  free  surface  reflections.  Note  that  this  converted  group  moves  rapidly  out  of 
the  synthetic  window  and  that  it  evolves  a  great  deal  ot  complexity  in  the  wavenumber  integration 
synthetics.  These  mean  that  it  will  not  strongly  affect  a  discriminant  focused  on  the  onset  of  Pg 
and  that  it  is  undesirably  sensitive  to  the  details  of  the  free  surface  gradient.  However,  it  may 
constitute  the  desirable  channel  for  scattered  S  energy  discussed  earlier.  At  any  rate,  the  fourth 
group  noted  in  Figure  10  is  a  second  head  wave  dominated  by  2Pn.  This  arrival  is  part  of  the 
same  generalized  ray  as  2PmP.  so  it  is  not  surprising  that  if  the  reflection  is  strong  so  is  the  head 
wave.  The  free  surface  reflection  plays  an  important  role  in  shaping  2Pn  just  as  it  does  for  2PmP. 
We  computed  similar  synthetics  for  a  crust  model  without  the  free  surface  gradient  and  found  that 
the  effect  of  the  gradient  is  to  dramatically  increase  the  significance  of  2PmP  and  all  of  the  free 
surface  reflections.  Without  the  gradient,  much  more  upgoing  eneigy  converts  to  SV  as  should  be 
expected. 

The  next  two  figures  present  those  records  which  we  believe  make  a  particularly  strong  case 
for  interpretation  of  the  Pg  observations  in  terms  of  the  four  basic  generalized  ray  groups  discussed 
in  the  previous  section  (PmP,  2PmP,  Converted  Group  and  2Pn).  Figure  I  I  shows  four  records  at 
distances  of  about  200  or  300  km.  Note  in  Figure  10  that  this  is  where  the  generalized  ray 
interpretation  is  most  valid  For  each  station,  we  show  the  observation  on  top.  a  complete 
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wavenumber  integration  synthetic  and  finally  a  generalized  ray  Green’s  function  beneath  to  aid  in 
interpretation.  The  two  records  from  MNV  comprise  a  strong  case  that  most  of  the  Pg  arrival  is 
comprised  of  an  interference  of  the  PmP  and  2PmP  group.  The  relative  amplitudes  and  exact 
details  differ  between  observation  and  prediction,  but  these  arc  two  second  waves  and  this  level 
of  fit  is  typical  of  that  generally  obtained  in  this  period  lange.  There  is  definite  support  for  the 
converted  ray  group,  but  its  role  in  determining  the  waveshape  is  not  as  strong  as  for  the  other 
two  groups.  The  records  at  K.NB  and  LAC  provide  further  suppoit  that  Pg  in  this  distance  range 
is  predominantly  an  interference  of  two  main  pulses  The  synthetics  again  do  not  predict  the  exact 
ratio  of  PmP,2PmP  amplitudes,  but  they  are  reasonably  close 

Figure  12  presents  a  similar  interpretation  of  four  records  at  close  to  400  km  range.  There 
is  continuing  evidence  of  the  importance  of  PmP  and  2PmP  in  determining  the  Pg  wave  shape. 
The  converted  group  is  outside  the  synthetic  window  The  primary  difference  is  that  the  2Pn 
group  has  emerged  in  front  of  PmP  and  caused  a  clear  pulse  it  is  important  to  note  the  inclusion 
of  the  PmS,  pPmS  and  pSmP  head  waves  within  the  2Pn  group.  The  reason  is  that  the  2Pn  wave 
itself  is  a  negative  pulse  as  shown  in  the  generalized  ray  Green’s  function.  In  all  cases  except  for 
the  El.K  record  the  first  motion  of  the  2Pn  group  is  clearly  positive.  The  wavenumber  synthetics 
predict  a  weak  positive  arrival.  We  interpret  the  strong  positive  arrivals  in  the  data  as  being 
associated  with  the  PmS  group,  though  the  observed  arrivals  are  definitely  sharper  and  higher 
amplitude  than  predicted.  Our  model  generally  fails  to  predict  head  wave  behavior  in  this  regard 
The  tit  of  the  observations  to  the  synthetics  after  the  onset  of  2Pn  is,  on  the  other  hand,  most 
acceptable  and  substantiates  our  model  tor  the  fine  structure  of  Pnl. 

The  adequacy  of  the  very  simple  crustal  structuie.  Model  3.  lor  modeling  the  completely 
trapped  (past  critical  angle;  response  though  not  the  head  waves  in  our  data  base  is  simple  to 
explain.  The  model  has  no  structure;  no  velocity  increase  at  all  beneath  the  crust-mantle  transition. 
The  completely  trapped  energy  has  no  penetration  beneath  the  Moho  so  it  is  insensitive  'o  this 
detail.  The  head  waves  Jo  penetrate,  and  m  particular,  would  be  sensitive  to  a  positive  gradient 
below  the  transition.  The  Moho  is  a  world-wide  feature,  so  that  :n  any  region  the  development 
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Figure  12.  Interpretation  of  regional  in  terms  of  the  major  ray  groups. 
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of  the  Pg  through  the  PmP-2l’mP  interference  could  be  anticipated,  and  because  the  eneigy  is 
completely  trapped,  the  waveform  would  be  insensitive  to  lateral  variations  in  lid  structure.  This 
has  important  implications  in  solving  the  event  discrimination  problem. 

Two  of  the  most  general  physical  bases  with  which  to  develop  a  discriminant  between 
explosions  and  earthquakes  are  the  event  depth  and  the  level  of  shear  radiation  We  have  found 
that  the  beginning  of  the  Pg  arrival  is  relatively  deterministic  and,  because  of  the  free  surface 
giadient,  lelatively  sensitive  to  the  free  surface  phases.  Because  of  the  increased  number  ol 
sigm: .cant  ray  paths,  the  number  ol'  free  surface  contributions  increases  along  with  the  complexity 
of  t1  .  .  .  rface  contribution.  Clearly  the  free  surface  interaction  is  very  dependent  on  source 

depth,  it  is  also  sensitive  to  the  level  of  shear  radiation  and  associated  S  to  P  conversion.  A 
clear  path  to  follow  in  attempting  to  develop  a  discriminant  would  be  to  locus  analysis  on  the  Pg 
onset  in  the  period  range  wo  have  been  studying  here.  To  establish  the  level  of  differences  to  be 
expected  between  explosions  and  earthquakes,  we  computed  Green’s  functions  for  a  point  double 
couple  buried  at  7  km  depth  in  our  Mode!  3.  The  p.cdictcd  waveforms  for  a  iarige  of  300  km  aie 
compared  to  the  synthetic  for  a  shallow  explosion  in  Figure  13  The  source  for  the  earthquake 
was  assumed  to  be  a  triangle  with  I  4  sec  rise  and  fall  time.  This  is  appropriate  for  a  magnitude 
5  b  earthquake  which  is  consistent  with  the  event  size  used  in  computing  the  explosion  source 
The  explosion  synthetic  is  shown  on  the  top  followed  by  synthetics  for  a  vertical  strike  slip,  a 
vertical  dip  slip  and  a  45°  dip  slip  fault.  The  differences  in  both  the  Pg  and  P„  wavofoinis  aie 
apparent. 

There  are  a  variety  of  approaches  possible  in  the  development  of  a  discriminant  from  the 
type  04  waveform  information  shown  in  the  figure.  In  the  case  of  Pn  we  used  a  simple  correlation 
norm  as  discussed  in  the  previous  section.  We  developed  a  discriminant  by  simply  determining 
the  average  waveform  expected  tor  an  explosion  and  correlating  that  waveform  with  Pn  waveforms 
from  a  mixed  data  base.  The  explosion-explosion  norms  were  on  average  about  30'Ni  higher  than 
the  explosion-earthquake.  The  correlation  norms  foi  the  earthquake  waveforms  with  the  explosion 
waveform  are  shown  next  to  the  traces  in  Figure  13.  The  correlation  of  the  explosion  with  itself 
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is  of  course  1.0,  but  it  is  reduced  by  at  least  40%  for  each  of  the  other  Green’s  (unctions.  From 
a  theoretical  standpoint,  discriminant  performance  comparable  to  the  Pn  discriminant  might  be 
anticipated.  Waveform  norms  were  computed  for  all  possible  source  mechanisms  by  stepping 
through  fault  orientation  angles  on  a  10"  grid.  The  average  correlation  decrease  of  the 
explosion-earthquake  was  45%  with  a  minimum  28%  For  noisy  signals  in  a  laterally  varying  earth, 
the  performance  of  a  correlation  discriminant  would  undoubtedly  be  worse  than  these  calculations 
indicate,  but  additional  investigation  seems  warranted. 

To  complete  this  preliminary  inquiry  into  the  potential  ot  the  Pg  waveform  time  domain 
discriminant,  we  have  computed  on  a  theoretical  basis  the  expected  trajectory  of  the  average  norm 
in  a  standard  discrimination  plot  (see  Figure  14).  To  compute  the  variation  of  a  typical  earthquake 
source  with  size  we  combined  relationships  of  seismic  moment  to  \1L  and  mb  (Taylor  el  a  I . ,  1988) 
with  those  for  source  duration  with  ^etsmiu  moment  for  western  North  America  (Somerville  et  al., 
I987j  to  obtain 

Logtdurutionj  =  5.77mb  -  2.7"8 

This  relation  has  the  same  information  requirements  as  the  Mueller- Murphy  source  which  we  are 
using  for  our  explosion  calculations  It  onlv  requires  an  mb  value  in  combination  with  a  depth 
scaling  relation  and  a  magnitude  yield  curve.  Figure  14  shows  the  trajectory  of  the  expected 
correlation  ot  an  explosion  Pg  waveform  with  an  earthquake  Pg  waveform  The  observed  P„ 
correlations  are  shown  lor  reterence  I  he  separation  of  the  observed  explosion  population  from 
the  theoretical  earthquake  population  is  lar  from  complete.  However,  there  is  clearlv  sufficient 
impetus  to  continue  the  investigations  ol  the  Pg  discriminant  further  and  every  reason  to  believe 
the  Pg  waveform  discriminant  will  be  as  good  as  any  developed  to  date.  We  again  emphasize  the 
physically  based  nature  (depth  phases)  of  the  Pn  and  Pg  discriminants  and  the  hope  this  offers  for 
diem  to  be  widely  transportable 
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Figure  14.  The  discrimination  performance  of  the  theoretical  cross  correlation 
of  Pg  expl .  to  earthquake  (crosses)  is  compared  to  t fie  observed  performance 
of  the  Pn  discriminant. 


tnfrnri»rfinn:  Understanding  the  detailed  structure  ot  the  seismic  discontinuities  of  the  Earth's 
upper  mantle  is  of  fundamental  importance  in  the  geosciences.  This  structure  and  its  regional  and 
lateral  variations  play  an  important  role  in  geodynamics  as  well  as  providing  constraints  for 
compositional  models  of  the  Earth  The  upper  mantle  P-wave  velocity  discontinuities  have  been 
studied  in  detail  for  many  different  regions  (e  g..  U.S.,  Canada.  Europe,  and  west  Asia)  However, 
until  recently,  data  recorded  in  central  Asia  have  not  been  readily  available,  preventing  the  Asian 
plateau  from  being  an  area  of  prolific  study. 

In  establishing  the  USA -USSR  Joint  Seismic  Program,  whereby  an  agreement  between  IRIS 
and  the  Soviet  Academy  of  Sciences  has  been  set  forth,  previously  unavailable  data  can  now  be 
obtained.  Decades  of  analog  records  from  the  Soviet  national  network  are  now  available,  and  are 
similar  in  many  respects  to  the  Woild  Wide  Seismographic  Station  Network  (Given,  Helmbergei 
and  Zhao,  1991.)  Many  records  from  this  database  are  already  digitized  and  easily  obtained  from 
the  Center  for  Seismic  Studies  (CSS/.  These  data  greatly  increase  the  number  of  event  to-station 
paths  in  central  Asia,  as  well  as  span  a  distance  range  bracketing  the  uppei  mantle  triplications, 
enabling  a  detailed  regional  study  of  the  upper  mantle  beneath  central  Asia. 

In  an  effort  to  study  the  upper  mantle  P-wave  velocity  profile  we  analyze  P-wave  motions 
at  distances  where  arrivals  from  the  upper  mantle  triplications  are  clearly  evident  The  data  used 
are  digitized  analog  vertical  component  short  period  records  of  USSR  underground  nuclear 
explusiuns.  A  comparative  approach  is  taken  wherein  USSR  data  are  correlated  with  predictions 
made  from  upper  mantle  models  presented  in  previous  studies.  These  models  were  derived  for 
distinctly  different  regions  elsew  here  on  the  globe,  and  are  used  here  in  an  effort  to  gain  a  qualitative 
understanding  of  the  lateral  variations  in  the  upper  mantle  of  the  different  locales  of  central  Asia. 
Region  of  Study.  A  map  of  the  study  area  is  presented  in  Figure  15.  The  eleven  stations  from 
which  we  have  short  period  analog  data  are  presented  as  triangles.  This  array  spans  distances  from 
around  5  to  40  degrees,  sampling  many  dilfer-nt  regions,  as  well  as  different  tectonic  settings. 
While  the  lithosphere  in  central  Asia  is  comprised  of  the  Russian  Platform  to  the  northwest,  and 


Figure  15.  Map  of  Asia.  1 1  stations  and  3  source  regions  used  Raypath  coverage  is  good, 
adequately  covering  central  Asia. 


the  West  Siberian  Platform  to  the  north,  the  regions  to  the  south  and  southeast  are  comprised  of 
complex  fold  systems.  The  question  of  whether  or  not  these  fold  systems  are  accompanied  by 
tectonic  rather  than  shield  seismic  velocities  in  the  lithosphere  is  important,  and  will  be  investigated 
in  this  study.  Also,  lying  between  the  West  Siberian  Platform  (plains)  to  the  north  and  the  Northeast 
Siberian  fold  system  (mountain  ranges)  to  the  east,  is  the  Central  Siberian  Platform  (plateau)  to 
the  northeast,  where  the  possibility  of  some  degree  of  a  tectonic  lithosphere,  or  transition  zone 
between  tectonic  and  shield  regions,  has  nr.f  yet  been  addressed  Three  stations  are  located  in  the 
southeast  (BOD,  TL  Y  and  TUP)  which  permits  us  to  investigate  such  a  possible  transition  zone. 
With  two  distinctly  different  and  densely  populated  source  regions,  Novava  Zemlya  and 
Sermpalatinsk,  there  are  many  path  geometries  across  Asia  to  investigate  such  lateral  variations 
(  omparaiive  Modeling  Technique:  A  compa.  Uive  modeling  approach  was  used  in  this  study, 
whereby  data  were  correlated  with  predictions  made  from  models  previously  presented  in  the 
literature.  The  synthetics  were  generated  by  the  generalized  ray  method  (Helmberger,  1973a).  The 
generalized  ray  technique  was  used  because  of  its  speed  and  cost-effectiveness.  The  procedure 
was  to  first  generate  Green's  functions  for  the  specific  models  and  distances  and  then  convolve 
them  with  empirical  sources  to  obtain  the  predictions.  Empirical  sources  were  estimated  by  taking 
several  P-wave  waveforms  at  teleseismic  distances  where  the  triplication  arrivals  are  well  separated 
in  time,  taking  care  to  choose  records  with  good  signal  to  noise  ratio  as  well  as  minimal  P-wave 
coda,  keeping  only  the  first  arrival  for  convolving  with  synthetics  (see  Figure  16).  Records  chosen 
were  typically  beyond  26  degrees  or  so,  a  distance  sufficiently  past  the  670  km  discontinuity 
triplication  crossover.  P-waves  with  different  source  time  functions  were  chosen  so  that  comparing 
synthetics  with  the  data  could  be  easily  accomplished  by  choosing  the  empirical  sources  having 
source  time  function  characteristics  matching  those  of  the  data  to  be  modeled. 

The  Data  Set:  Our  data  set  is  a  collection  of  479  digitized  short  period  vertical  component  records 
from  II  stations  of  the  Soviet  national  network,  recording  I  10  underground  nuclear  explosions. 
The  number  of  events  recorded  at  eacn  station  is  presented  in  Figure  17a.  Some  stations  are  very 
densely  populated,  while  others  are  sparse  The  distance  distribution  of  the  events  recorded  for 


Figure  16. Synthetic  Seismogram  Technique.  An  empirical  source  time  function  is 
derived  bv  taperinu  the  P  wave  of  a  record  where  no  triplications  arrivals  exist  in  the 
time  window  of  interest.  This  is  then  convolved  with  the  time  derivative  ot  a  general¬ 
ized  ray  generated  Green’s  function  to  obtain  the  model  prediction. 
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Short.  Period  Analog  Records 


Figure  1 7  la)  Population  o.  events  recorded  at  each  station.  The  average  is  about  4 
recordings  per  station,  (b)  The  number  of  recordings  plotted  with  respect  to  distance 
Tnplication  distances  are  well  populated  for  this  data  set. 
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this  network  is  shown  in  Figure  17b.  As  this  figure  illustrates,  the  distances  where  the  different 
triplication  arrivals  are  routinely  recorded  from  around  15  to  27  degrees  are  well  populated.  Nearly 
all  of  the  110  events  are  from  the  two  main  Soviet  test  sites,  Semipalatinsk  and  Novayn  Zemlya, 
with  the  exception  of  several  off  site  Peaceful  Nuclear  Explosions  (PNE’s.)  The  Novaya  Zemlya 
events  are  from  the  northern  Matochkin  site.  The  Semipalatinsk  events  are  from  3  specific  sites: 
Degelen  Mountain,  Kovnstan  and  Shagan  River.  The  Semipalatinsk  events  are  plotted  in  Figure 
IS  using  the  USGS  PDE  locations.  The  data  from  each  of  the  three  Semipalatinsk  sites  were  anlyzed 
separately.  The  CK.M-3  instruments  recording  these  events  ate  similar  in  nature  to  the  World  Wide 
Seismographic  Station  Network  (WWSSN)  short  period  instruments,  peaking  at  around  1  second. 
This  data  set  is  easily  obtained  in  digital  form  from  CSS. 

The  controlled  nature  of  the  source-receiver  geometries  along  with  the  abundance  of  the 
short  period  data  makes  it  ideal  to  stack  the  data.  Records  were  binned  with  respect  to  source 
regions  and  stations,  and  a  separate  stack  was  made  for  each  station  at  each  source  site,  provided 
sufficient  data  existed  for  any  given  path.  Separate  stacks  were  made  for  the  three  specific 
Semipalatinsk  sites  shown  in  Figure  18.  Thus,  for  any  given  station,  up  to  four  stacks  may  have 
been  made:  one  for  each  of  the  Semipalatinsk  sites,  and  one  for  the  northern  Novaya  Zemlya  site. 
The  data  for  specific  paths  contain  similar  waveforms,  and  after  stacking,  the  triplication  signals 
become  enhanced  while  coda  signals  decrease  in  amplitude.  The  distance  window  for  any  site-receiver 
path  usually  varies  no  more  than  about  30  km,  so  the  stacks  were  made  by  lining  up  the  first 
arrivals  of  the  waveforms.  All  amplitudes  were  normalized  to  unity  before  stacking.  An  example 
for  the  well  populated  station  OBN  is  given  in  Figute  19.  The  distance  window',  using  the  USGS 
PDE  locations,  spans  only  about  7  km.  The  resulting  stacked  traces  of  the  short  period  data  set 
are  used  to  model  upper  mantle  structure.  Records  were  stacked  only  for  distances  where  triplication 
arrivals  are  clearly  seen  above  the  noise  level.  Records  were  not  included  in  a  stack  if  they  had 
portions  that  were  zero  in  amplitude  presumably  due  to  undigitizable  sections  of  the  original 
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Figure  18.  Location  of  events  from  three  Semipoiaiinsk  sites.  The  locations  are  from 
the  USGS  PDE.  The  data  from  each  site  were  stacked  separately. 
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Figure  19.  Stacking  Technique.  For  a  specific  source  region  and  station,  the  first  ar¬ 
rivals  are  lined  up,  all  amplitudes  are  normalized  to  the  peak  amplitude,  then  the  traces 
are  summed.  The  example  given,  for  station  OBN  recording  Shagan  River  events,  the 
variation  in  distance  spans  about  10  km.  The  resulting  sum  exhibits  enhanced  triplica¬ 
tion  arrivals  and  suppressed  coda  arrivals. 
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seismograms.  The  total  number  of  records  fitting  the  above  criteria  is  183.  For  the  binning  of 
records  chosen  above,  20  stacks  have  been  made,  with  the  average  number  of  records  per  stack 
being  about  9  (Table  1). 

Models  for  Comparative  Study  Synthetic  seismograms  were  generated  for  the  distances  in  Table 
1  for  four  different  models.  These  synthetics  were  then  compared  to  the  data  stacks  in  an  effort 
to  obtain  a  first  order  idea  of  the  type  ot  model  necessary  to  explain  the  data  for  each  path.  The 
models  chosen  represent  different  types  of  lithospheric  structure  (Figure  20),  such  as  a  thick  shield 
model  (S25,  LeFevre  and  Helmberger.  1990,  Canadian  shield),  a  model  with  no  low  velocity  zone 
(KCA,  King  and  Calcagmle,  1976,  west  Russia),  and  two  models  with  a  low  velocity  zone  in  the 
lithosphere,  K.8  (Given  and  Helmberger,  1980,  northwestern  Eurasia)  and  T7  (Burdick  and 
Helmberger,  1978,  western  U.S.),  where  T7  is  a  tectonic  model.  Although  these  four  models  were 
derived  for  different  regions  of  the  globe,  the  comparison  of  theii  predictions  to  the  Soviet  data 
is  important  to  assess  the  possible  nature  of  the  underlying  upper  mantle.  The  most  important 
differences  in  these  models  are  in  die  top  200  km.  Second  order  differences  are  the  depths  of  the 
"400”  and  “670“  km  discontinuities,  gradients  above  and  below  the  discontinuities,  as  well  as  the 
percent  velocity  increase  of  the  discontinuities.  For  the  rest  of  this  paper  we  loosely  use  the  numbers 
"400"  and  "670”  to  represent  the  two  major  upper  mantle  discontinuities,  keeping  in  mind  that  their 
exact  depth  varies  from  model  to  model,  as  well  as  in  the  Earth. 

Absolute  TraveL Times:  The  arrival  times  of  the  first  arrivals  of  the  large  short  period  data  set 
were  usually  very  easy  to  measure  due  to  the  first  arrival’s  impulsive  nature.  Records  having  low 
signal  to  noise  ratio  were  not  used  in  our  travel  time  analysis.  One  of  the  uncertainties  in  the 
travel  times  is  the  use  ot  the  USGS  PDE  origin  times  (unfortunately,  the  origin  times  are  not 
available  from  the  USSR.)  Also,  the  records  came  from  CSS  with  a  file  of  instrument  time 
corrections  that  range  trom  0  seconds  up  to  around  60  seconds.  Up  to  three  time  corrections  per 
day  are  listed  for  each  station  for  the  days  that  events  occurred.  In  many  cases,  the  clock  error 
was  not  constant  over  the  period  of  a  day,  and  changed  by  as  much  as  a  second  and  a  hall  This 


Table  1:  Stacked  Traces  List 


STACK 

NO. 

STA 

NAME 

DELTA 

(deg) 

# 

REGS 

SOURCE 

REGION 

SITE 

NAME 

1 

ARl) 

13.43 

3 

Semipalaunsk 

Degelen  Mountain 

2 

ARU 

13.75 

4 

Semipalatinsk 

Shagan  River 

3 

TLY 

15.72 

14 

Semipalaunsk 

Shagan  River 

4 

ARU 

17.09 

9 

Novaya  Zemlya 

Matochkin  Shar 

5 

OBN 

19.76 

7 

Novaya  Zemlya 

Matochkin  Shar 

6 

NR1 

20.14 

4 

Semipalatinsk 

Koynstan 

7 

NVS 

21.86 

8 

Novaya  Zemlya 

Matochkin  Shar 

8 

BOD 

22.01 

4 

Semipalatinsk 

Shagan  River 

9 

BOD 

22.49 

3 

Semipalatinsk 

Degelen  Mountain 

10 

NVS 

22.87 

4 

W.  Kazakh 

PNE 

11 

APA 

24.36 

2 

Semipalatinsk 

Degelen  Mountain 

12 

APA 

24.60 

"> 

Semipalatinsk 

Shagan  River 

13 

OBN 

25.22 

5 

Semipalatinsk 

Koynstan 

14 

TUP 

25.28 

5 

Semipalatinsk 

Shagan  River 

15 

OBN 

25.56 

17 

Semipalatinsk 

Degelen  Mountain 

16 

OBN 

25.96 

39 

Semipalatinsk 

Shagan  River 

17 

BOD 

27.44 

9 

Novaya  Zemlya 

Matochkin  Shar 

18 

UZH 

28.31 

5 

Novaya  Zemlya 

Matochkin  Shar 

19 

TLY 

29.88 

5 

Novaya  Zemlya 

Matochkin  Shar 

20 

UZH 

36.20 

31 

Semipalatinsk 

Shagan  River 

35 


Depth  (km) 


Mode 
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Canada 
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UP  Body 
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W  Russia 
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NW  F.urasia 
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Body  waves 
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Burdick  <£  Helmberget,  1978 

UP  Body 
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Figured.  Velocity  Models.  Four  models  were  used  :n  this  study  fur  comparison  with 
the  data  stacks.  The  most  outstanding  differences  in  the  models  are  in  the  uppennost 
200  km  of  the  mantle.  S28  has  a  fast  thick  lid  and  a  gradual  low  velocity  zone,  17 
and  K8  have  pronounced  low  velocity  zones,  and  KCA  has  no  low  velocity 


zone. 


may  be  another  source  of  travel  time  error.  For  cases  where  the  clock  error  changed  within  a  day, 
we  interpolated  the  given  time  corrections  for  that  day  to  the  P-wave  arrival  time  for  that  day's 
event,  assuming  the  clock  error  changed  in  a  linear  fashion 

First  arrival  travel  times  are  shown  in  Figure  2)  with  the  predictions  of  model  KCA. 
Anomalies  are  clearly  visible  on  this  figure,  especially  at  15.’’  degrees  (station  TLY)  where  the 
first  arrival  is  as  delayed  as  the  discontinuity  reflection  This  nr>ay  be  a  slow  tectonic  path  from 
the  Semipalatmsk  region. 

With  the  uncertainty  in  travel  times  due  to  origin  time,  along  with  the  data  being  grouped 
at  '60  discrete  distances  as  seen  in  Figure  21,  a  ~  investigation  is  not  possible.  Furthermore, 
because  of  these  uncertainties,  we  choose  not  to  use  absolute  travel  times  from  this  data  set  as  a 
discriminant  for  choosing  between  the  different  models. 

Waveform  Comparisons  As  Figure  15  indicates,  our  study  area  has  good  coverage  of  the  ~entral 
Asia  area.  However,  the  source- receiver  geometry  is  such  that  no  given  path  from  a  source  region 
has  more  than  one  station  on  it.  Without  a  profile  of  stations  for  a  specific  path,  obtaining  a 
unique  solution  model  is  not  possible.  While  deriving  a  definitive  model  is  not  practical  in  this 
situation,  comparisons  with  model  predictions  can  still  be  made  to  assess  agreement  with  different 
classes  of  models  and  assess  lateral  variation  from  path  to  path. 

Predictions  for  the  four  models  K8,  KCA,  S25  and  FT  are  compared  to  the  data  stacks  in 
Figure  22.  Nine  panels  are  presented,  in  each  of  which  the  top  trace  is  the  data  stack,  and  the 
four  underlying  traces  are  model  predictions.  The  time  window  shown  is  30  sec,  and  the  prominent 
arrivals  are  a  function  of  distance.  An  example  of  the  relative  timing  order  of  the  triplication 
arrivals  is  seen  in  the  travel  time  curve  for  KCA  in  Figure  21.  L’p  to  around  21°  the  first  arrival 
is  the  direct  P  from  above  tl  e  400  (PAB),  and  tbe  second  arrival  is  the  combination  of  the  reflection 
off  the  top  of  the  400  km  discontinuity  (PBC)  and  the  P  wave  from  below  the  400  (PCD).  *n  some 
cases,  the  reflection  off  of  the  top  of  the  670  (PBE)  along  with  the  direct  arrival  from  below  the 
670  (Pp,f)  are  evident  as  a  third  arrival.  For  these  distances,  the  goodness  of  fit  criteria  are  how 
well  wavefoims,  amplitude  ratios  and  differential  travel  times  in  the  synthetics  agree  with  the 
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Figure  21.  First  Arrival  Travel  Times.  First  arrival  picks  of  the  short  period  data 
made  from  impulsive  waveforms  only  are  shown  here  with  the  prediction  of  model 
KCA.  The  letters  A  through  F  represent  the  names  of  the  different  triplication 
branches. 
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f  igure  Z2(a)  through  (i)  present  (i;ita  stacks  nop  trace)  and  the  predictions  of  the  four 
models  (Ixittom  four  traces.)  Comparisons  are  described  in  the  text 


data.  This  distance  window  corresponds  to  Figure  22  (a)  through  (d).  From  right  before  the 
crossover  distance  of  the  PAB  and  PCd.  to  a  little  after  the  crossover  of  the  PCd  and  the  arrival 
from  below  the  670  (PEF),  the  differential  time  between  the  triplication  arrivals  is  often  too  small 
to  measure.  This  corresponds  to  Figure  22  (e)  through  (h),  where  at  least  two  arrivals  appear  as 
more  of  a  wave  packet  than  distinct  arrivals.  The  shape  ot  these  wave  packets  (data  and  synthetics) 
are  very  sensitive  to  the  amplitude  ratios  of  the  different  arrivals,  and  thus  can  be  a  good  diagnostic 
tool  to  distinguish  between  different  possible  models.  Beyond  the  crossover  distance  ot  P(:I)  and 
PEF  the  arrivals  separate  in  time  so  that  differential  travel  time  information  again  supplements  the 
amplitude  ratios. 

The  degree  of  fit  of  the  model  predictions  to  the  data  stacks  varies  from  good  to  poor.  In 
some  cases  differential  travel  time  information  is  well  predicted,  but  amplitude  behavior  is  not. 
Also,  the  short  period  data  is  undoubtedly  affected  by  site  sti  ucture  reverberation,  which  we  make 
no  effort  to  account  for  here  othei  than  attempting  to  suppress  this  affect  in  the  stacking  process 
Nonetheless,  the  triplication  arrivals  are  usually  strong  and  clear  in  this  data  set.  Data  from  many 
of  tne  paths  to  the  west  of  the  test  sites  in  Figure  1  .5  are  better  modeled  in  the  first  order  sense 
by  a  shield  or  shield-type  model.  Data  to  TLY  to  the  east  from  Semipalatinsk,  however,  is  better 
modeled  by  the  tectonic  western  US  model  T7  Figure  22  (b)  shows  the  data  stack  and  predictions 
tor  station  TLY  at  15. 7°.  Only  model  T7  predicts  the  anomalously  small  differential  travel  time 
between  PAB  and  Pod  due  to  its  slow  lid.  To  a  slightly  more  northerly  azimuth  is  station  BOD. 
Figure  22  has  two  central  Asian  stacks  for  BOD.  With  about  a  0.5  degree  difference  in  distance. 
Figure  72  ( g j  is  for  the  Shagan  River  site,  while  Figure  22  (h)  is  for  the  Degelen  Mountain  site. 
For  Figure  22  tg),  the  prediction  from  K.8  is  the  best.  Increasing  in  distance  0.5°  (Figure  22  ( h ) ) 
yields  a  record  that  is  poorly  modeled  by  all  lour  models  It  >s  possible  that  the  path  to  BOD  from 
Semipalatinsk  may  coincide  with  a  transition  region  or  boundary  between  the  tectonic  region  to 
the  south  and  a  tnote  shield- like  lithusphere  to  the  north,  producing  complex  records. 


Station  OUN  recorded  many  events  from  the  three  Sermi'alatinsk  sites.  The  three  resulting 
013N  data  stacks  and  model  predictions  are  presented  in  F  igure  2.3  For  this  distance,  the  first 
arrival  is  PEp  and  the  second  arrival  is  a  combination  of  Pod  and  Pde-  Around  8  or  10  sec  after 
the  first  arrival,  the  back  branch  of  the  400  triplication  consisting  of  PAB  and  Pgc  can  be  seen  in 
the  synthetics.  A  remarkable  feature  in  the  data  in  Figure  23  is  the  rapid  diminishing  with  distance 
of  the  second  arrival  (Pcd  and  Pde)  This  feature  is  not  predicted  in  the  synthetics.  Given  and 
Helmberger  (1980)  presented  data  in  this  same  distance  window  from  Semipalatinsk  to  an  azimuth 
to  the  southwest  recorded  at  TAB  (Tabriz,  Iran  )  For  that  azimuth,  the  changing  amplitude  ratio 
of  the  first  two  arrivals  is  well  predicted  by  K8  (Figure  24)  "The  OBN  data  appears  anomalous  in 
that  the  second  arrival  decays  so  rapidly  in  such  a  small  distance  window  of  around  0.8°.  This 
feature  of  diminishing  PCD  and  PDE  relative  to  P|.;f  in  a  small  distance  window'  (i.e.,  move  the  tip 
ot  the  triplication  to  a  smaller  distance)  may  be  more  accurately  modeled  by  increasing  the  velocity 
gradient  between  the  400  and  670  discontinuities.  Increasing  the  gradient  causes  down-going  rays 
in  this  depth  range  to  turn  up  more  sharply.  Thus  the  deepest  possible  ray  propagating  above  the 
discontinuity,  i.e  ,  the  arrival  at  the  tip  of  the  D  branch,  is  bent  up  to  a  smaller  distance. 

Station  TUP  has  an  azimuth  from  Semipalatinsk  intermediate  to  that  of  stations  TLY  and 
BOD  for  that  same  source  region.  For  this  path,  the  distance  to  TUP  is  25.28°.  This  distance  is 
in  the  window  of  that  presented  in  Figure  23  to  station  OBN,  though  the  waveforms  are  quite 
different.  The  source  time  function  of  the  TUP  record  is  longer  period,  which  may  partly  obscure 
any  differentiation  between  the  two  triplication  phases  that  are  arriving  just  a  few  seconds  apart. 
Also,  if  these  waves  propagate  along  a  transition  zone  from  a  tectonic  lithosphere  to  that  of  a 
shield,  complexities  should  be  expected.  Figure  25  (a)  displays  the  stack  for  TUP  and  the  synthetics. 
The  prediction  from  model  K.8  gives  the  closest  fit,  though  it  is  still  poor.  The  prediction  from 
the  tectonic  model  T7  is  not  as  good  as  that  of  K8.  It  is  possible  that  the  transition  from  tectonic 
to  shield  is  abrupt  because  the  T7  model  is  the  only  model  that  adequately  predicts  the  record  from 
the  neighboring  path  to  the  south  at  station  TLY  (Figure  22  (b)).  For  paths  to  TLY,  TUP  and  BOD 
from  Semipalatinsk,  the  models  K.8  and  T7  produce  motions  closest  to  our  observations.  This 
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Figure  23.  OBN  from  Scm.palatmsk  The  three  columns  are  separate  slacks  ...  i"c  ' 
sues  in  Figure  4  and  model  |'i c.licrions  lor  llial  speeme  .nslanec  For  I  rest  . 
the  rellection  off  the  670  arrives  around  7  or  .«  seconds  al.er  the  lira  arrival,  and 
arrival  of  die  400  arrives  around  7  or  so  seconds  after  the  first  arrival  (For  tins  mv 
lanee  range,  only  S75  predicts  a  visible  -«KJ  arrival.  and  "  -  'x-ry  small  Tl~s 
receiver  pair  shows  no  strong  sign  ot  an  amval  m>m  die  ■  )  .)  ■  “ 
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Figure  24.  Three  Semipalatinsk  events  recorded  at  TAB  (Tabriz,  Iran)  are  shown  here 
with  K8  synthetics  (from  Given  and  Helmberger,  1980.)  The  amplitude  ratio  of  the 
first  two  arrivals  is  well  predicted  by  K8  tor  this  southwest  azimuth  from  Semipala- 
tinsk. 
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Figure  25.  (a)  through  (c)  present  data  stacks  (top  trace)  and  the  predictions  of  the  tour 
models  (bottom  four  traces.)  Comparisons  are  described  in  the  text. 
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implies  that  over  a  lateral  distance  of  around  500  km  or  so,  the  low  velocity  zone  might  drastically 
change  from  one  like  that  in  T7  to  one  like  that  in  K8.  Finer  modeling  with  broadband  data  is 
necessary  to  resolve  such  issues. 

For  larger  distances  (Figure  25  (b)  and  (c ) )  the  predictions  widely  vary  for  the  four  models 
presented  here.  The  observed  amplitude  ratio  of  the  first  two  arrivals  in  each  of  (b)  and  (c)  is 
not  predicted  by  the  synthetics.  Also,  models  KCA  and  S25  predict  a  large  late  arrival  that  has 
reflected  off  of  the  400  km  discontinuity.  This  is  not  seen  in  the  data.  If  the  large  second  arrival 
m  the  data  is  PCD  with  PDE,  then  S25  does  better  than  the  other  models  in  predicting  the  differential 
travel  time,  though  the  amplitude  ratio  prediction  is  very  poor  Modifications  in  the  fine  details 
of  the  gradients  above  and  below  the  discontinuities  are  necessary  to  model  these  two  paths. 
Discussion:  The  four  models  chosen  for  this  comparative  study  represent  widely  varying  upper 
mantle  rheologies.  All  four  models  were  derived  from  body  wave  studies.  KCA  is  a  logical  choice 
as  a  comparison  model  because  it  was  derived  from  short  period  data  in  the  western  Russia  region. 
K  8  was  chosen  because  it  has  a  shield-like  lithosphere  as  well  as  being  derived  for  parts  of  central 
Asia  that  coincide  with  out  study  area.  K8  was  constructed  using  KCA  as  a  starting  model, 
utilizing  both  short  and  long  period  data  from  Soviet  explosions.  A  classic  shield  model  is  S25, 
derived  for  the  Canadian  shield  region  using  long  period  earthquake  data.  A  first  order  justification 
in  choosing  S25,  a  model  that  was  derived  for  a  completely  different  locale,  is  that  central  Asia  is 
thought  to  be  shield-like,  with  the  East  European  Craton  to  the  west  and  the  Siberian  Craton  to 
the  east.  This  choice  has  proven  useful,  as  S25  is  a  best-fitting  model  for  several  of  the  paths  in 
our  study  area.  To  assess  regions  where  tectonics  are  thought  to  play  an  important  role,  model 
T7  was  utilized.  T7  was  derived  for  the  tectonic  western  US  fron  long  period  body  waves  of 
earthquakes  and  explosions.  Other  models  in  the  literature  may  have  done  equally  well  for  our 
criteria  of  choosing  four  models  with  fundamentally  differing  lithospheres.  We  note  again  that 
our  choice  of  models  is  to  emphasize  regional  and  lateral  variations  in  the  upper  mantle  beneath 
Asia,  and  not  to  make  an  argument  for  specific  models  for  the  specific  locales. 
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We  have  not  used  absolute  travel  time  information  in  our  modeling  procedure.  Due  to 
uncertainties  in  hypocenter  information,  we  are  unable  to  analyze  travel  times  at  well  populated 
stations  with  array  data  techniques.  However,  we  can  still  make  qualitative  comparisons  between 
the  different  stations.  For  example,  data  recorded  at  TLY  from  Semipalatinsk  are  anomalously  slow 
by  up  to  8  seconds  when  compared  to  KCA  predictions  (Figure  21,  15.7°).  Absolute  travel  times 
for  the  direct  P  waves  at  this  range  for  the  various  models  show  considerable  differences.  The 
largest  difference  is  between  S25  (pure  shield)  and  GCA  (pure  tectonic,  Waick  ( 1 084 ))  is  roughly 
13  seconds.  The  travel  time  for  GCa  at  15.7“  is  over  10  seconds  slower  than  that  of  KCA,  which 
brackets  the  TLY  times  in  Figure  7.  T7  predicts  an  arrival  time  near  3  seconds  earlier  than  the 
TLY  data,  though  modifications  in  the  direction  of  a  more  tectonic  model  makes  matching  these 
observations  possible. 

The  13  source-receiver  paths  that  are  most  relevant  to  upper  mantle  triplication  distances 
are  presented  in  Figure  26  along  with  the  best-fitting  model  for  that  path.  For  some  paths,  two 
models  are  listed  because  predictions  from  either  model  may  be  considered  as  a  better  fit  than  the 
others  The  model  name  is  in  parentheses  if  the  best-fitting  model  for  that  path  poorlv  predicts 
observed  motions.  If  no  model  came  close  to  predicting  observations,  a  question  mark  is  shown. 
A  comparative  study  has  been  presented  in  this  report,  so  that  the  four  models  in  Figure  26  are 
primarily  meant  to  give  a  first  order  view  of  the  upper  mantle  regimes  beneath  Asia,  as  well  as 
possible  lateral  variations  from  iocale  to  locale.  All  of  the  models  used  in  this  study  have  a  LVZ, 
with  the  exception  of  KCA  Only  one  path  favored  KCA  (station  OBN  from  Novaya  Zemlya), 
with  predictions  from  S25  being  equally  satisfactory.  A  first  order  qualitative  conclusion  from  this 
comparative  approach  is  that  from  the  four  models  chosen  for  comparison  in  this  study,  the  north 
and  the  west  regions  of  our  study  aiea  favor  shield  (S25)  or  shield-type  (K8)  models,  while  the 
area  to  the  southeast  is  probably  a  transition  zone  from  a  shield-type  mantle  to  a  more  tectonic 
upper  mantle  (T7).  A  more  detailed  modeling  approach  is  necessary  for  each  path  to  make  more 
quantitative  inferences  about  the  structure  of  the  underlying  mantle. 


Figure  26.  Best  fitting  models  for  this  data  set  are  written  on  the  wave-paths.  If  the 
best  fitting  model  or  models  do  a  poor  job  at  best  tn  predicting  observations,  the  model 
name  is  put  in  parentheses 
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Several  paths  in  this  data  set  produce  motions  that  are  poorly  predicted  by  the  above  models. 
For  example,  the  paths  from  Semipalatinsk  westward  to  stations  ARU  and  OBN  are  close  in  azimuth 
and  are  poorly  modeled.  For  ARU,  the  differential  travel  time  and  amplitude  ratio  of  PAg  and 
PCD  >n  data  are  not  predicted  in  any  of  the  synthetics  (Figure  22a).  For  OBN,  rapid  decay  of 
the  second  arrival  in  a  very  shcrt  distance  window  is  seen  in  the  data  and  not  the  synthetics  (Figure 
23). 

In  the  case  of  station  OBN,  the  diminishing  second  arrival  at  25.96“  can  be  attributed  to  an 
early  ending  of  the  “D"  tip  of  the  670  km  triplication  To  produce  such  a  feature  in  the  synthetics, 
changes  must  be  made  to  the  model  in  the  transition  zone  Mode)  changes  made  far  above  the  400 
km  discontinuity  tend  not  to  diminish  the  amplitude  of  the  second  arrival  enough  to  match  the 
data.  Support  for  this  statement  is  seen  in  Fdgure  23  in  that  four  very  different  lithospheric 
structures  fail  to  produce  a  diminished  second  arrival  One  way  to  end  the  670  km  triplication  at 
a  smaller  distance  would  be  to  include  a  small  zone  of  an  increased  velocity  gradient  on  the  top 
side  of  the  670  discontinuity  This  turns  rays  up  more  sharply,  thus  decreasing  the  distance  where 
the  "D"  tip  occurs  An  example  of  such  a  scenario  is  presented  in  Figure  27.  The  left  column  ot 
synthetics  is  from  K8,  whereas  the  middle  column  is  from  K.8.I,  where  K8  has  been  modified  by 
a  small  zone  of  increased  gradient  right  above  the  670  discontinuity.  The  prediction  of  K.8.1  at 
26. 0‘  adequately  suppresses  the  second  arrival.  It  has  been  recently  proposed  that  the  depth  of  the 
670  km  discontinuity  may  \arv  by  as  much  as  20  km  Note  also  that  the  depth  of  this  discontinuity 
for  different  regions  of  the  globe  has  been  placed  anywhere  from  around  640  km  to  690  km  in 
the  literature  A  670  km  discontinuity  occurring  at  a  more  shallow  depth  would  also  give  rise  to 
the  "D"  tip  of  the  triplication  ending  sooner.  An  example  ol  the  discontinuity  raised  by  17  km 
along  with  a  slightlv  increased  overlying  gradient  is  presented  in  the  third  column  of  Figure  77 
This  model  K.8.2  also  predicts  the  correct  amplitude  behavior  at  26  0°.  These  types  of  foiward 
modeling  calculations  can  be  carried  out  at  each  station.  Although  we  note  here  that  the 
non-uniqueness  due  to  the  source-receiver  geometry  prevents  the  validation  of  such  results. 


Figure  27  Predictions  are  plotted  from  K.X  and  two  models  with  increased  transition 
rone  gradients,  where  one  has  a  more  shallow  670  discontinuity  by  17  km.  The  ampli¬ 
tude  ratio  of  the  first  two  arrivals  for  K8.1  and  KS.2  is  opposite  of  that  for  KX.  This 
amplitude  behavior  is  seen  in  the  data  at  OBN 


49 


Recently,  there  has  been  renewed  debate  over  the  existence  of  a  520  km  discontinuity. 
Synthetics  were  made  for  a  model  D25  t  Jones  and  Helmberger  199l\  a  m  idel  where  S25  was 
modified  by  adding  a  1  percent  velocity  increase  at  500  km  depth,  and  to  preserve  travel  time 
the  increase  in  velocity  at  670  km  was  decreased.  In  comparing  D25  synthetics  with  the  data,  a 
discontinuity  of  this  nature  is  not  supported  by  the  short  period  data  stacks  in  this  study. 

With  broadband  data  from  the  Soviet  Union  becoming  more  available,  future  work  tor  this 
region  will  include  a  detailed  forward  modeling  effort.  By  using  other  stations  in  the  Soviet  national 
network  that  were  not  .ncluded  in  the  CSS  data  set,  along  with  selected  earthquakes,  a  better  view 
of  the  Asian  upper  mantle  may  be  constructed.  For  the  source -receiver  geometry  inherent  in  this 
data  set,  wave  paths  are  typically  isolated  such  that  no  multi-station  profiles  exist  for  any  path. 
For  this  reason  we  have  not  explored  solution  models  for  any  of  the  data  due  to  non  uniqueness. 


COMPOSITION  Qf  REGIONAL  Pn  and  Sn  WAVEFORMS  RECORDED 

IN  THF.  I  .S.  and  IJ.S.S.R 


Introduction:  With  the  methodologies  available  to  compute  high-frequency  full-wave  responses 
both  for  vertically  and  laterally  inhomogeneous  crustal  media,  the  broadband  modeling  of  regional 
displacement  seismograms  has  become  a  topic  in  many  recent  studies  for  retrieving  the  source  and 
propagation  path  characteristics  (Barker.  1991,  Burdick  e!  a!..  1991;  Helmberger  el  al..  1991;  Zhao 
and  Helmberger,  1991).  The  method  of  Cagmard  de-Hoop  is  generally  used  to  understand  the 
composition  of  various  seismic  phases  that  the  source  process  and  the  propagation  effect  together 
make  up  at  a  receiving  station.  This  method  requires  many  generalized  rays  to  be  tracked  between 
the  source  and  the  receiver  As  the  crustal  medium  becomes  complicated,  the  method  can  quickly 
become  quite  cumbersome  with  the  process  of  just  tracking  the  rays.  For  this  reason  alone,  the 
method  developed  based  on  the  frequency  -  wavenumber  integration/  reflectivity  approaches  is  highly 
useful.  This  later  method  is  computation  intensive,  but  with  the  advent  of  fast  computers  and 
vectorized  machines  this  is  no  longer  a  constraint  on  the  methodology. 

At  high  frequencies  the  effects  of  scattering  in  the  crust  become  so  intense  that  only  statistical 
properties  of  wave  forms  are  meaningful.  However,  one  exception  to  this  generality  involves  the 
long-period  Pn  and  shear-coupled  P  waves  at  regional  (l°-15°)  ranges  (Helmberger,  1 9^2 ,  1973). 
In  a  recent  study,  Saikia  and  Burdick  (1990)  showed  that  the  short-period  Pn|  waves  (period  as 
short  as  2  s)  are  also  stable.  They  studied  many  observations  from  the  Nevada  Test  Site  (NTS) 
explosions  recorded  at  regional  distances  of  200  to  420  km  and  modeled  the  Pn|  waveforms  using 
a  deterministic  crustal  waveguide  The  sources  of  these  waveforms  were  shallow.  Also,  the  sources 
were  predominantly  isotropic,  and  the  portion  of  Pn|  waves  which  was  included  in  the  Pn  and  Pg 
waves  had  a  duration  of  about  30  s  and  was  dominated  by  compressional  waves.  To  understand 
the  observed  data,  Saikia  and  Burdick  (1990)  employed  the  following  strategy.  They  used  the 
frequency- wavenumber  algorithm  iCODH;  FILON_aS,  written  by  Chandan  K.  Saikia, 
Woodward-Clyde  Consultants)  to  compute  the  explosion  generated  Pni  waves  for  several  canonical 
crustal  models  and  selected  a  crustal  model  based  on  the  agreement  between  the  data  and  the 


synthetic  seismograms.  The  crustal  model  was  then  utilized  to  understand  the  composition  of  the 
Pg  wave  group  which  was  constituted  of  phases  like  PmP,  pPmP,  2PmP,  PmS,  pPmS,  PmPSmP, 
PmSPmP  etc.  Pg  is  a  constituent  wave  group  whose  frequency  content  is  widely  used  to  discriminate 
events.  In  this  study,  we  have  taken  a  similar  strategy  to  investigate  the  broadband  composition 
of  Pni  and  Sn|  seismograms  recorded  in  the  North  American  continent  at  iegional  distance  from 
double-couple  sources  and  of  the  P  and  S  waveforms  that  are  recorded  within  the  Soviet  Union 
We  shall  mainly  focus  on  identifying  the  rays  important  to  model  the  regional  waves  within  the  S 
wave  window. 

Data:  For  the  U.S.  study,  we  used  a  set  of  three-componen'  broadband  seismograms  recorded  at 
Harvard  (HRV)  station  at  a  distance  of  640  km  from  the  Saguenay  earthquake  of  November  25, 
1988  (Figure  28a).  These  seismograms  were  recoroed  on  a  Streckeisen  seismometer.  We  selected 
these  seismograms  because  many  features  recorded  on  the  seismograms  were  successfully  modeled 
by  Zhao  and  Helmberger  (1991).  Beginning  with  their  crustal  model,  we  have  directed  our  study 
towards  the  modeling  the  high-frequency  details  observed  in  the  Pn|  waves  and  the  composition 
of  waves  identified  as  Sn  and  sSn  by  Zhao  and  Helmberget  (1991)  using  a  multiple  source  model 
A  similar  study  was  directed  towards  the  modeling  of  tne  broadband  seismograms  recorded  within 
the  Susie t  Union.  We  selected  a  set  of  three-component  seismograms  recorded  at  GARNI  from  an 
earthquake  which  originated  at  a  distance  of  200  km  at  an  azimuth  of  290u  on  May  4,  1989  (38.73°N 
and  "8.5°\V,  Figure  28b).  Unlike  for  North  American  earthquakes,  the  waveforms  of  very  few 
Soviet  Union  earthquakes  have  been  modeled  Thus,  it  is  necessary  to  develop  a  starting  crustal 
model  even  to  obtain  a  first-order  agreement  between  data  and  synthetic. 

Modeling  of  HR  V  Seismograms  from  the  Saguenay  Lai  thqual.e:  Figure  29  shows  the  broadband 
displacements  recorded  at  Harvard  station.  lo  investigate  the  influence  of  crustal  structure  on  tin- 
various  significant  phases  of  the  Pni  window,  we  have  started  with  the  crustal  model  shown  in 
Figure  >0  t.  \  Jutted  lines,  this  model  extends  from  the  surface  to  a  half  space  at  a  depth  of  55 
krn.  The  mator  velocity  discontinuity  is  at  a  depth  of  55  km  where  the  P  velocity  pimps  from 
0  "  I  km  s  to  8  i  km.-  5  and  the  S  velocity  jumps  I  torn  5  .52  kin.  s  to  4  V  k  m.  s.  /. Iiao  and  I  lelmbci  per 


Figure  28.  Geographical  location  of  (a)  HRV  (Harvard)  station  (solid  triangle)  and  November  25, 
1988  Saguenay  earthquake  (solid  star)  and  (b)  GAR  (Garm)  station  (solid  triangle)  and  May  4,  1989 
USSR  earthquake  (solid  star). 
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Broadband  Displacement  recorded  at  Harvard 
Station  from  1988-11-25,  Saguenay  Earthquake 
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Figure  29.  Broadband  three-component  displacement  seismograms  as  recorded  by  Harvard  station 
trom  the  1938.  November  25  Saguenay  earthquake.  The  original  seismograms  were  integrated 


ENA  Regional  Velocity  models 


Figuie3Q.  Preliminary  regional  crustal  model  developed  for  tl.e  region  between  the  Garm  station 
and  the  source  of  the  1989  May  4,  USSR  earthquake 


(1991)  used  a  reflectivity  code  (Mallick  and  Frazer,  1988)  to  compute  the  medium  response  and 
used  an  elastic  crustal  structure  to  model  the  data.  They  used  a  y,, (shear- wave  quality  t actor)  ot 
6200  and  stated  that  a  lower  value  of  is  not  required  to  match  the  recorded  wave  form,  although 
the  conventional  wisdom  is  that  for  eastern  North  America  On  is  of  the  order  of  300  (Hwang  and 
Mitchell,  1987).  Based  on  this  published  information,  we  started  to  look  for  certain  phases  within 
the  Pn;  regime  for  which  the  agreement  between  the  data  and  synthetic  can  be  improved  and  in 
the  process  to  learn  more  about  the  legional  waveguide. 

The  phases  marked  as  Sn  and  sSn  show  the  greatest  misfit  between  the  data  and  the  synthetics 
computed  by  Zhao  and  Helmberger  (1991)  (see  their  Figure  16).  The  synthetic  seismograms  are 
definitely  of  lower  frequency.  So  our  initial  attempt  was  to  understand  what  part  of  the  crustal 
waveguide  would  be  most  critical  in  development  of  these  waveforms.  In  the  present  calculation, 
we  used  the  frequency-wavenumbei  integration  method  and  set  the  nyquist  frequency  at  10  hz. 
We  computed  theoretical  seismograms  for  eight  fundamental  faults  and  used  a  focal  mechanism 
with  a  dip  65°,  a  rake  of  78c  and  a  strike  of  323°  to  predict  the  vertical,  radial  and  tangential 
component  seismograms.  These  synthetics  were  used  to  compute  both  the  point  and  multiple  source 
seismograms  and  the  corresponding  vertical  component  seismograms  are  shown  in  Figure  31.  The 
source  model  contained  three  sub-sources,  with  seismic  moments  of  1. 55x10*4,  I.45xl024  and 
1 .9i x  I O24  dyne-cm  respectively.  The  second  source  was  delayed  by  0.65  s  and  the  third  source 
by  1.45  s  from  the  first  source  to  account  for  the  propagation  of  the  rupture  front.  The  first 
source  was  represented  with  a  source  time  function  defined  by  a  trapezoid  of  0  4s  rise  time,  0.05s 
of  follow-on  time  and  0.25s  of  healing  time.  Similarly,  the  second  and  third  sources  were  convolved 
with  the  trapezoids  of  (0.2s,  0.15s,  0.15s)  and  (0  Is.  0.3s.  0.2s),  respectively.  W'e  also  show  the 
synthetic  seismograms  generated  by  Zhao  and  Helmberger  (1991)  in  Figure  32  using  a  nyquist 
frequency  of  4  hz  so  that  a  direct  comparison  can  be  made  with  those  shown  in  Figure  31.  The 
frequency  content  in  the  Pnl  waves  uf  these  seismograms  is  not  as  rich  as  those  Pnl  waves  shown 
in  Figure  31  In  these  seismograms  (Figure  31),  the  high  frequencies  are  the  result  of  derived 
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Point-Source  Vs  Multiple-Source  Synthesis  of  the 
Harvard  Record  for  Various  Crustal  Models 
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Figure  31.  Comparison  between  two  sets  ot  displacement  seis.nogratns  synthesized  using  point  and 
multiple  sources,  (a)  Point-source  displacements  -  the  upper  seismogram  is  computed  using  the 
mouei  response  of  Zhao  and  Helmberger  (1991)  and  the  botto.;  seismogram  is  computed  using  the 
model  response  of  the  present  crustal  model,  and  ib.)  multiple -source  displacements  lor  the  two 
crustal  models. 
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Figure  32.  Comparison  between  data  and  synthetic  displacements  with  a  nvquist  of  4  hz  The 
seismograms  for  625  km  was  computed  using  a  different  velocity  crustal  model  (Figure  taken  from 
Zhao  and  Helmberger,  1991) 


In  Figure  33,  we  compare  the  vertical  and  radial  component  showing  just  the  Pnl  portion  of 
the  seismograms  computed  using  the  parameters  of  multiple  sources.  The  high-frequency  signals 
are  adequately  predicted  with  respect  to  those  observed  on  the  recorded  data.  The  seismograms 
computed  using  the  response  up  to  4  Hz  were  essentially  identical  to  these  seismograms. 

Rai  Analysis  of  Pnl  Seismograms:  In  this  section,  we  discuss  our  investigation  of  the  constituent 
phases  of  the  recorded  Pnl  seismogram  at  Harvard  station.  The  basic  idea  is  to  investigate  the 
interaction  of  individual  ray  groups  in  creating  the  total  seismogram.  We  computed  generalized 
ray  seismograms  using  the  source  process  of  the  Saguenay  earthquake  for  several  groups  of  generalized 
tays  In  Figure  34,  we  display  vertical-component  seismograms  of  these  ray  groups.  The  top  six 
seismograms  are  normalized  to  their  maximum  amplitude.  All  the  PmP  and  SmS  rays  were  allowed 
to  reflect  from  each  interface  beneath  the  crust-mantle  boundary  including  the  reflection  from 
the  Moho  discontinuity.  The  total  response  of  these  PmP  and  SmS  rays  is  plotted  in  the  first 
seismogram.  The  geometric  arrivals  are  indicated  by  PmP  and  SmS  respectively.  The  Sn  arrival 
is  small  and  is  preceded  by  a  refracted  phase  SP.  This  refracted  phase  had  developed  due  to  a 
critical  incidence  of  an  S  wave  on  an  interface  permitting  the  converted  P  phase  to  travel  along 
the  interface.  This  is  a  strong  phase  as  distinctly  observed  on  the  record  (see  Figure  33)  and  was 
incorrectly  identified  by  Zhao  and  Helmberger  (199!)  as  the  Sn  wave.  The  seismogram  in  the 
second  row  is  for  the  sPmP,  a  ray  which  has  departed  from  the  source  as  a  S  wave  and  then 
converted  to  P  mode  at  the  free  surface  The  amplitude  of  this  ray  is  small  The  next  seismogram 
is  tor  sSmS.  Both  the  geometric  and  head  waves  are  strong  for  this  ray  group  and  contribute 
significantly  to  the  total  seismogram.  The  next  two  seismograms  are  for  the  SmSSmS  and  sSmSSmS 
ray  groups.  Both  the  ray  groups  have  significant  contributions.  The  sixth  seismogram  is  for  a  ray 
group  identified  as  SmS'SmS.  The  rays  included  in  this  group  leave  the  source  downward  and 
reflect  from  each  interface.  The  reflections  ate  turned  back  into  the  lower  crust  again  at  the  Moho 
discontinuity  before  they  are  reflected  back  to  the  receiver.  The  contributions  from  these  rays  do 
ofer  a  significant  contribution  to  the  evolution  of  the  Sn  wave  group.  The  seismogram  "Total"  is 
the  result  ot  direct  sum  of  the  upper  six  seismograms.  Having  obtained  a  good  agreement  between 
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Figure  }3.  Comparison  between  data  and  synthetic  displacem  mts  computed  with  model  response 
i-P  to  a  nyquist  of  10  Hz  Note  the  de'.elopmeni  of  higti  t  requeue ies  and  agreement  between  phases 
marked  by  the  z-rows. 
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Figure  34.  Understanding  of  the  waveform  recorded  at  Harvard  station  using  the  ray  decomposition 
technique.  The  top  six  seismograms  are  for  the  individual  tay  groups.  The  seismogram  labelled 
"Total"  is  the  total  response  of  all  the  responses  of  upper  six  seismograms  and  the  comparison  with 
the  F-K  seismogram  shown  below  suggests  a  good  agreement  between  the  two  seismograms 
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the  data  and  the  synthetics,  we  only  plotted  the  multiple-source  frequency  wavenumber  seismogram 
computed  using  the  frequency-wavenumber  method  beneath  the  total  response  for  a  direct 
comparison.  This  comparison  produced  good  agreement  among  the  dominant  features  within  the 
so  called  "Sn  waves". 

Thus,  we  have  extended  our  previous  study  (Saikia  and  Burdick,  1990)  on  the  deciphering 
of  the  ray  composition  of  Pg  waves  from  explosion  sources  to  earthquake  sources.  As  in  the  above 
study,  we  found  that  the  waveforms  within  the  Sn  group  can  oe  studied  in  time  domain  in  terms 
of  a  basic  few  rays,  namely  the  SmS,  sSmS,  SmSSmS,  sSmSSmS  and  SmS'SmS  rays.  Since  these 
phases  leave  the  source  as  S  waves,  they  are  not  excited  by  the  explosion  source.  Therefore,  the 
only  phases  that  may  arrive  within  the  Sn  widow  from  a  pure  isotropic  source  are  the  P  waves  that 
are  converted  to  S  waves. 

Mode'ing  Qf  GARNI  Seismograms  from  Mav  4»  1989  USSR  Earthquake:  In  this  study  we  have 
used  a  set  of  three-component  seismograms  recorded  at  Garm  station  from  an  earthquake  of  May 
4,  1989  (latitude:  39  436°N  and  longitude:  75  35°F,  h  =  35  km,  ISC).  The  station  is  located  at  a 
distance  of  200  km  from  the  source  Figure  35  shows  the  recorded  displacement  processed  from 
the  broadband  velocity  seismograms.  A  high-pass  fitter  was  applied  to  remove  the  long-period 
effects  The  crustal  structure  encountered  by  the  wavefield  along  its  propagation  path  is  complex 
which  is  reflected  in  the  waveform.  To  begin  to  understand  the  waveforms,  it  was  necessary  to 
develop  a  crustal  structure. 

Our  strategy  for  developing  the  crustal  model  was  to  begin  with  the  tangential  component 
seismogram  because  of  the  simplicity  of  the  observed  displacement.  This  component  contains  only 
three  distinct  individual  arrivals  as  marked  bv  the  arrows.  We  used  generalized  ray  theory  to 
synthesize  this  component  for  various  crustai  models  and  several  source  depths.  The  best  prediction 
was  obtained  for  a  source  depth  of  25  km  for  the  preliminary  crustal  model  shown  in  Figure  36. 
We  used  the  following  fucal  mechanism,  dip-75^  s I i p=  -  1 3 5°  and  a  strike=32°.  The  structure 
contains  two  major  discontinuities  representing  the  Conrad  and  the  Moho.  We  succeeded  in  modeling 
the  tangential  displacement  using  only  three  arrivals,  the  direct  SH  and  two  reflections  from  the 


Broodband  Displacement  recorded  at  GARM  Station 
from  1989-05-04  Earthquake  in  Soviet  Union;  AZ  =  292 
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Tigure  35.  Broadband  three-component  displacement  seismograms  as  recorded  by  Garm  station 
from  the  USSR  earthquake  of  May  4,  1989  The  original  seismograms  were  integrated. 
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USSR  Regional  Velocity  Model 


Figure  3d.  tasiern  North  American  regional  velocity  profiles.  The  models  shown  by  the  solid 
lines  are  developed  m  this  stud>  The  model  Jiown  be  the  dotted  lines  were  published  b>  Zhao 
and  Helmbercer  ( l‘>(>  1 ). 


discontinuities.  The  frequency- wavenumber  seismograms  were  computed  with  these  parameters 
and  were  compared  with  the  recorded  data.  Figure  37  shows  the  comparison  between  the  dat3  and 
the  synthetic  seismograms  for  the  vertical  and  tangential  components.  We  are  successful  :n  producing 
a  good  agreement  between  the  data  and  the  synthetic  lor  the  tangential  motion.  We  marl,  the 
inu.viduai  arrivals  in  the  synthetic  and  show-  their  correspondence  with  the  data  by  the  thin  arrows 
I  he  vertical  component  show  agreement  in  the  arrivals  of  Pn.  PmP,  ScS  and  SmS  phases.  The 
signal  bracketed  within  the  window  of  the  vertical-component  synthetic  seismogram  has  a  similar 
character  in  the  frequency  content  to  the  signal  bracketed  within  the  data  window.  Figure  38 
shows  a  generalized  ray  seismogram  using  the  thiee  rays.  The  vertical  synthetic  seismogram  has 
a  strong  SmS  which  is  smeared  out  in  the  data  due  to  the  interaction  with  the  physically  more 
complicated  crust  in  the  region.  The  vertical  component  of  the  recorded  seismogram  is  also 
dominated  by  long-period  Pnl  signals  shown  by  the  solid  window.  These  waves  are  also  observed 
in  the  synthetic  seismogram  but  arriving  at  Ge.rm  with  a  fast  velocity.  We  also  investigated  the 
effect  of  a  possible  linear  velocity  gradient  near  the  free  surface  to  determine  if  suen  velocity 
distribution  would  account  for  the  mismatch  between  the  data  and  the  synthetics  withe,  picketed 
(he  window  shown  in  Figure  37  We  discretized  the  top  ten  kilomeicrs  of  the  crust  into  ten  layers 
of  equal  thickto  .nd  allowed  a  P-wave  velocity  increase  from  4.5  km. /sec  to  5.5  km  sec  from 
the  surface.  The  S-wave  velocity  within  each  layer  had  a  ratio  of  1.73  to  the  P-wave  velocity. 
This  seemed  to  be  a  particularly  reasonable  explanation  lor  the  small  complex  phases  between  the 
majoi  arrivals.  However,  the  synthetic  seismograms  computed  using  this  surface  gradient  did  not 
improve  the  fit  to  a  significant  degree. 

G'lnp.u  isoii  of  Regional  Pnl  Waves  from  die  Lii  and  LciMi  kJJiiidi  Models:  In  this  section,  we 
continue  to  investigate  the  regional  waveforms  tint  are  likely  to  be  predicted  by  the  crustal  models 
developed  for  North  America  and  Soviet  Union.  Since  the  recoid  at  Harvard  station  was  so 
successfully  modelled  and  since  it  was  at  a  lange  of  640  km,  wc  examined  the  response  of  the 
USSR  crust  model  at  this  range.  In  I  ict,  seismograms  at  such  distance*  are  just  becoming  available 


the  Soviet  Union. 


Figure  30  shows  the  comparison  between  me  synthetic  seismograms 


Preliminary  Modeling  of  Broadband  Displacement 
Recorded  at  GARM  station.  R  =  200  Km  I 
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Figure  38  A  vertical  component  seismogram  computed  using  direct  P  and  S,  PcP  and  bcS.  an 
I’mP  and  SmS  phases.  These  arrivals  can  distinctly  be  obser'ed  on  the  recorded  seismogram. 
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Comparison  of  Displacement  Seismograms  at  640  Km 
for  the  U.S.  and  USSR  Crustal  Models 
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Figure-  39.  Comparison  between  lire  synthetic  seismograms  computed  31  640  km  for  the  USSR 
crustal  model  and  the  US  crustal  model  Both  vertical  and  radial  components  are  shown. 


computed  for  the  USSR  and  US  crustal  models  for  both  the  verti  ents  assuming 

the  same  focal  mechanism  at  the  same  azimuth  The  seismograms  ^  at  the  respective 

depths  of  the  earthquakes  and  the  two  depths  are  similar.  We  also  used  the  same  source  function. 
The  USSR  crustal  model  predicts  a  stronger  PmP  relative  to  the  Pn.  The  Sn/SP  and  sSn  arrivals 
predicted  in  the  seismograms  by  the  US  crustal  model  (marked  by  the  arrows)  seem  to  exhibit  a 
correspondence  to  the  long-period  signals  predicted  in  the  response  of  the  USSR  crustal  model. 
The  difference  in  the  amplitude  ratios  is  caused  bv  the  differences  in  the  near-surface  velocities 
of  the  two  crustal  models.  For  the  North  American  crustal  model,  the  crustal  velocities  have  a 
gradient  near  the  surface.  The  rays  arrive  at  the  receiver  more  steeply  compared  to  the  rays  for 
the  USSR  crustal  model,  thus  partitioning  the  energy  in  a  significantly  different  ratio  to  the  vertical 
and  radial  component. 

We  further  investigated  the  composition  of  the  Pnt  waves  predicted  by  the  USSR  crustal 
model  at  640  km  in  terms  of  generalized  rays.  We  found  that  the  S-wave  reflections  from  the 
Mulio  and  Conrad  discontinuities  are  strong  as  shown  in  the  top  two  seismograms  of  Figure  40. 
The  Conrad  reflection,  ScS,  arrives  immediately  following  the  Moho  reflection  SmS.  The  phase 
shown  by  an  arrow  on  the  SmS  seismogram  is  arriving  at  the  arrival  time  of  Sn  phase,  but  its 
waveshape  is  more  complicated  than  is  expected  from  a  classical  Sn  phase.  Among  the  other 
phases  that  contribute  most  significantly  to  the  total  Pn)  seismogram  within  the  S  window  are  the 
sSmS,  SmSSmS  and  ScSScS  phases.  The  bottom  two  seismograms  plotted  in  Figure  40  allow  us  to 
compare  the  agreement  between  the  generalized  ray  (Total)  and  frequency  wavenumber  (F-K.) 
seismograms.  The  agreement  is  poor  following  the  PmP  arrival.  A  strong  long-period  signal  does 
propagate  to  the  receiver  in  the  frequency  wavenumber  seismogram.  This  must  be  a  total  effect 
of  many  generalized  rays.  This  effect  was  also  observed  on  the  recorded  seismogram  3t  Garm 
station  even  at  a  distance  of  200  km  from  the  source. 

Conclusions:  Based  on  these  investigations,  it  seems  feasible  to  develop  time-domain  discriminants 
at  different  nuclear  test  sites  which  rely  on  the  stable  features  that  are  observed  in  the  recorded 
waveforms.  The  most  stable  phases  are  observed  in  the  explosion  generated  Pn|  waveforms  for 
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Figure  J0.  Understanding  of  the  waveform  computed  at  640  Km  from  the  USSR  crustal  model 
using  vne  ray  decomposition  technique  The  top  six  seismograms  are  for  the  individual  rav  groups 


periods  as  short  as  2  s  (Burdick  et  at..  1988;  Saikia  and  Burdick,  1990).  In  this  studs,  we  have 


extended  our  analysis  approach  to  regional  broadband  seismograms  from  earthquake  sources  The 
short -period  Pn|  waves  have  a  functional  dependence  on  the  crustal  waveguide.  They  can  be 
deterministically  modeled  using  average  flat- layered  ctustal  structures  and  using  some  selected 
generalized  rays.  By  modeling  the  broadband  displacement  at  Harvard  station,  we  found  that  the 
structure  across  the  crust-mantle  transition  zone  and  within  the  mantle  can  profoundly  affect  the 
Irequency  content  of  the  phases  like  Sn  and  sSn.  The  source  multiplicity  of  an  earthquake-  can  aiso 
cieate  added  complexity  in  the  frequency  content  of  these  phases  We  found  that  the  Pn|  seismograms 
near  the  S-wave  arrival  can  adequately  be  modeled  using  the  ray  responses  of  the  following  phases: 
SmS,  sSmS,  SmSSmS,  sSmSSmS  and  SmS'SmS. 

For  the  Soviet  Union,  the  most  important  requirement  for  understanding  recorded  seismograms 
is  the  crustal  model.  The  structure  within  the  Soviet  Union  is  heterogeneous  and  the  development 
of  reliable  crustal  models  is  on-going  (Gurrola  and  Minster.  1991 ).  In  this  study,  we  have  developed 
a  crustal  structure  by  modeling  the  recorded  seismograms  at  Garm  station  from  an  earthquake  at 
a  distance  of  200  km  (Az=292°)  which  consisted  of  Conrad  and  Moho  discontinuities.  In  addition, 
a  slight  gradient  is  allowed  for  the  uppei -mantle  structure.  The  rav  analysis  indicated  that  the 
most  important  generalized  rays  for  the  composition  of  the  Pn]  waves  are  the  following  phases: 
SmS,  ScS,  sSmS,  sScS,  SmSSmS  and  sSmSSmS. 
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